General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



'v/eA 


I 



u 

O' 

< 

LU 

CO 


GT 


> 


CO 

o: 

LU 



DEPARTMENT CF PHYSICS 

SCHOOL OF SCIENCES AND HEALTH PROFESSIONS 

OLD DOMINION UNIVERSITY 

NORFOLK, VIRGINIA 


TECHNICAL REPORT NO. PTR-83-5 

A STUDY OF PHYSICAL PROCESSES 
FOR SPACE RADIATION PROTECTION 

By 

G.S. Khandelwal 
and 

Y.J. Xu 



Progress Report 

For the period ending Decenber 1983 


O 

2 


Prepared for tne 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, Virginia 




Under 

Research Grant NCCI-42 

Dr. Lawrence W. Townsend, Technical Nbnitor 
, Space Technologies Branch 

P 

# (NASA-CIt- 173 I 83) A STUDY OF PHYSICAL 
j PROCESSES FOR SPACE RADIATION PROTECTION 
Progress Report, period ending Dec. 1903 
(Old Dominion Univ., Norfolk, Va.) 108 p 
HC A 06/ttF A 0 1 CSCL 20H 






January 1984 


DEPARTMENT OF PHYSICS 

SCHOOL OF SCIENCES AND HEALTH PROFESSIONS 

OLD DOMINION UNIVERSITY 

NORFOLK, VIRGINIA 


TECHNICAL REPORT NO. PTR-83-5 


A STUDY OF PHYSICAL PROCESSES 
FOR SPACE RADIATION PROTECTION 


By 

G. S. Khande^wal 
and 

Y.J. Xu 


Progress Report 

For the period ending December 1983 


Prepared for the 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, Virginia 


Under 

Research Grant NCC)-42 

Dr. Lawrence W. Townsend, Technical Monitor 

Space Technologies Branch 




Submitted by the 

Old Dominion University Research Foundation 

P.0. Box 6369 

Norfolk, Virginia 23508 


January 1984 


A STUDY OF PHYSICAL PROCESSES FOR SPACE RADIATION PROTECTION 

By 

G.S. Khandelwal 1 and Y. J. Xu 2 

The performance on the above project, during this period, can best be 
described by the following three items (two research papers and one 
dissertation), which are described in detail in the papers attached with 
this report. 


1. Intermediate Energy Proton Stopping Power for Hydrogen Molecules 
and Monoatomic Helium Gas, Physics Lett . A (in press). 


2. Low-energy proton stopping power of N 2 * 0 2 , and water vapor 
(Revised for publication in the Physical Review.) 


3. A theoretical model for calculating molecular stopping power, 

Y. J. Xu's dissertation, to be submitted to the faculty in partial 
fulfillment of the requirement of the Ph.D. degree in Applied 
Physics. 


T Eminent Professor, Department of Physics, Old Dominion University, Norfolk, 
Virginia 23508. 


2 Graduate Research Assistant, Department of Physics, Old Dominion 
University, Norfolk, Virginia 23508. 
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Intermediate Energy Proton Stopping Power 
for Hydrogen Molecules and Monoatomic Helium Gas 
Y. J. Xu and G. S. Khandelwal 
Department of Physics# Old Dominion University 
Norfolk# Virginia 23S08 
U.S.A. 
and 

J. W. Wilson 

NASA Langley Research Center 
Hampton# Virginia 23665 
U.S.A. 


ABSTRACT 

Stopping power in the intermediate energy region (100 keV to 1 MeV) 
has been investigated, based on the work of Lindhard and Winther, and on 
the local plasma model. The theory is applied to calculate stopping power 
of hydrogen molecules and helium gas for protons of energy ranging from 
100 keV to 2.5 MeV. Agreement with the experimental data is found to be 


within 10 percent 


Recent investigations suggest the usefulness of the local plasma 
model of Lindhard and Scharff (1) in calculating the mean excitation 
energies [2-5] of various systems. As is known, such calculations of mean 
excitation energy help in readily obtaining stopping power for high energy 
incident particles. On the other hand, the so called intermediate energy 
region (100 keV to 1 MeV in the context of this paper) is of considerable 
interest. Unfortunately, the theoretical calculations in this energy 
region are hindered, primarily due to the complexity of the problem. 
Bonderup's [6] theoretical calculations of stopping power for some 
elements, however, seem to be in good agreement with experimental results 
for Incident proton energy above 500 keV. 

In this letter, we present calculations of stopping power based on 
the work of Lindhard and Winther [7] , and on the local plasma model of 
Lindhard and Scharff. The theory is expected to cover the intermediate 
energy region while retaining its applicability for the high energy 
case. The method should be easily applicable for calculating stopping 
power of different materials composed of molecules such as H 2 and 
gases, water vapor, and many other chemical systems of interest. The 
method first establishes a stopping number function L, based on Lindhard 
and Winther' s formulas. Next, the Gordon-Kim model [8] for density 
function of a h 2 molecule is used to calculate its stopping power in the 
energy region from 100 keV to 2.5 MeV. Calculations are also performed 
for the monoatomic Helium gas. 
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0RIG1NAL RAGS 
OF POOR QUALITY 


The celebrated stopping power formula for an energetic charged 
particle of charge ze and velocity v , traversing matter of charge number Z 
is given by: 

dE 4Tiz 2 e‘ 4 ... 

“jr— NZL (1) 


dx mv 

where m is the mass of an electron, and N is the number of atoms per unit 
volume of the medium. 

The stopping number L of Eq. (1) has been a topic of considerable 

» 

study. Most interest, however, has usually been centered on the high 

2mv 2 


energy case for which L attains a simple formula, L - In 


In this 


situation, only one parameter I, the mean excitation energy of the medium, 

determines the stopping power. Lindhard and Winther have investigated the 

function L for a free electron gas in the regions of low and high energy 

incident charged particles. These authors [6, 7] give for the high energy 

1 


case, the expression for L to order — 2 as: 


L = In Y - 


< T > 

1 o 
— mv‘ 


(2) 


2 mv 2 4 2 

where, Y = , the classical plasma frequency w 2 = — , p the 

p pm 

electron density, and < T > , the average kinetic energy is given by 

< T > - mv F 2 ' 
where v^ is the Fermi velocity. 

For the low energy case, they give: 

.2 3 /4 3 /2 


Y z J/4 3 /t 

(f-> * '~yx) 


(3) 


where. 


with 


c 1 (x) 


2 (1 - j ) 2 


log 


2 = 


' * ix 2 


1 ~ ?* 2 

1 * 4y 2 


X = 


Tihv. 


The calculational procedure in this paper proceeds with the 

exploitation of Eqs. (2) and (3) with the following modifications in 

mind. First, as shown in a previous reference in context with the 

calculations of mean excitation energy, due to the Pine correction, one 

should replace til by yoo in the equations for L function where, y is given 
P P 


by [2, 9] 


' s 


» “ 1 + f t <’ + T? P 2 ' 

s 


(4) 


where \ * — ^ is the average distance between electrons in units of Bohr 

3 a o 

radius. The parameter (3 is determined by solving the following equation 

0.076 p 2 + - 0.916 = 0 

A. / 


The equation (2) for the L function warrants some discussion. First, 
one notes that the L function of Eqs. (2) and (3) are derived by Lindhard 
and Winther for a free electron system. Transition to an atomic system of 
the first term of Eq. (2), as studied widely, is accomplished under the so 
called local plasma model in which density p (r) is evaluated by using 
Quantum mechanical wave functions. Such a Quantum mechanical analogue of 
the second term of Eq. {%) would be of interest. In this context, a 
result first derived by Brown [10] would prove to be useful. Brown 
studied the K-shell asymptotic stopping power of an Hydrogenic system 
(with 2 K-electrons) for a fast projectile, taking the maximum momentum 
transfer to a free electron appropriately equal to 2mv. He obtained for 

the second term in Eq. (2) for an Hydrogenic system, the quantity 

1 < T > 

— — — ■ — - . Thus, m order to make a transition to an atomic system 

2 1 2 

2 mV 


consistant with the free electron model, we assume Brown's result and 
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accordingly replace Eq. (2) by: 


L = In Y - 


3 1 * 5 1 


(5) 


10 X Y 

At this stage, we combine the low and high energy L functions to 


determine the appropriate dependence of stopping power on energy. To do 
thi3, we used equations (3), (4) and (5) for our desired results. 

Bonderup, in his calculations, had combined Eqs. (2) and (3) and assumed a 
constant value of y equal to /2. Furthermore, we attempt to preserve 
continuity between (r, v) , the low energy stopping number function 
given by Eq. (3) and (r, v) , the high energy function given by Eq. 

(5), in the following way. Extensive numerical evaluations, using a 


computer program, of the functions and for various values of the 


variable Y, revealed that in most useful cases there were found two 


roots Y^ and Y^ of the equation L^= where Y^ < Y^« Furthermore, the 


slope M of the function for values of Y greater than Y^ was always 


small. On the other hand, the slope of for Y values less than Y^ was 


very steep. These observations, including the behavior of these functions 
(see the figure), led us to the following recipe to preserve the 
continuity consistant with the physics of the situation: 


L - L 


when 


MY < L 


L « Max [L , L 3 


when 


2 

MY > L, 


(6) 


where, and L 2 are given by Eq. (3) and (5) respectively. Eq. (6) can 
now be used to obtain stopping number values for a system, given the 
velocity of the projectile and the density p (r). For a molecule, Gordon- 
Kim give the density as: 

p <r) - p, <r> ♦ p 2 [r - R > (7) 


' f 


i •„ 


> - ♦ 

}i* * 4 -* , 


« * 
. ; * " 

if •*■!.*, 

ei . . v.*‘/ 
- . * * * 

* . 

k \ * • 


K* 

k * 




where, (r) is given by the square of the ground state wave function of 

the hydrogen atom. R 12 is the distance between two hydrogen atoms, which 

is known to be 1.4a Q . Knowing the density by Eq. (7), stopping 

power was obtained by integrating Eq. (6) over r, for protons of energy 

between 100 keV and 2.5 MeV, for both H 2 and monoatomic helium gas. For 

He atoms, Hartree-Fock wave functions [11] were employed. Table 1 lists 

the results of this paper, together with Anderson-Ziegler curve fitted 

results r 1 2 } , and three sets of experimental data [13, 14, 15]. Good 

agreement, within 10 percent, is found with the experimental data. 

Calculations of the equations established by Bonderup were also undertaken 

by extending them to molecular systems. Table 1 lists these values for 

the H molecule in the last column. These differ from our results in the 
2 

low energy region. Table 2 lists the same physical quantities for helium 
gas [14-18]. The same trend is observed as in the case of the H^olecule. 
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stopping cross sections . (in units of ev x lo” cm /atom) 

for a Heli.ua atom for various proton energies 



* Anderson and Ziegler [12) note that the data in reference [14] ■ rt ^ P -r 
constantly high by amount ranging fro* 5 percent to 10 percent. 





































FIGURE CAPTION 


The low and high energy stopping number functions as functions of the 
variable Y. A typical value equal to 0.1 of the quantity was chosen. 
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ABSTRACT 

A modified local plasma model, based on the works of Lindhard and 
Winther; Bethe, Brown, and Walske, is established. The Gordon-Kim model for 

molecular electron density is used to calculate stopping power of N 2 , O 2 , 
and water vapor for protons of energy ranging from 40 keV to 2.5 MeV, 
resulting in good agreement with experimental data. Deviations from Bragg's 
rule are evaluated and are discussed under the present theoretical model. 


2 


Recently, departures from Bragg's rule have been noticed in the theo- 
retical calculations of the mean excitation energies of various molecular 
i 2 

systems. * Analysis of the experimental data on energy loss of low energy 
a particles in gases also indicates deviations from Bragg's rule, 3 - 5 In 
this paper, the stopping power theory of Lindhard and Winther, 6 and the 
local plasma theory of Linoiiard and Scharff, 7 are used to perform calcula- 
tions in the low energy region. Modifications are introduced through a 
simplifying model which incorporates the effects of the shell corrections 
and of the screening of the projectile. The model is justified on the basis 
of fulfilling the more ambitious aim of obtaining the molecular stopping 
power. The Gordon-Kim 8 electron density model of molecular wave functions 
is utilized in the calculations. Such a model, as is known, allows a 
successful method of calculating chemical bond effects. Calculations done 
on N 2 , Og and water vapor are found to be in fair agreement with experi- 
ments. Furthermore, departures from the Bragg's rule are noticed for all 
these systems. 

The celebrated stopping power formula for an energetic charged particle 
of charge ze and velocity v, traversing matter of charge nunber Z is given 


by 


dE _ 4ir j^e 4 
dx mv 2 


NZL , 


( 1 ) 


vrtiere m is the mass of an electron, and N is the nunber of atoms per 
unit volume of the medium. 

The stopping nunber L of Eq. (1) has been a topic of considerable 
study. For instance, Lindhard and Winther have investigated the function L 


3 


for a free electron gas in the regions of low and high energy incident 
charged particles. These authors give for the high energy case, the expres 
sion for L to order as 


In Y - 


< T > 
j mv 1 2 3 


( 2 ) 


where, Y = , the classical plasma frequency u> p = J > p 

the electron density, and < T >, the average kinetic energy is given by 

< T > - (jjj) mvp 2 , 

where Vp is the Fermi velocity. 

For the low energy case, they give 


i * (5 ) 


2 3/4 


3/2 

Y 



where 


Mx) = 


1 

v 2 2 
(i -•*-) 

3 


log. 


1 + ! x 2 



(3) 



with 


Equation (2) for the L function warrants some discussion. First, one 
notes that the L function of Eqs. (2) and (3) are derived by Lindhard and 
Winther for free electron system. Transition to an atomic system of the 
first term of Eq. (2), as studied widely, is accomplished under the so- 


called local plasma model in vtfnch density p (r) is evaluated by using 
quantun mechanical wave functions. The local plasma model is equivalent to 

replacing the molecular dipole oscillator strengths by the corresponding 
classical plasma absorption spectrum. The adequacy of such a replacement 
was recently shown by Johnson and Inokuti 9 to be most accurate for evaluat- 
ing atomic quantities associated with stopping power inspite of differences 
between the plasma spectrim and the actual oscillator strength distribution. 
A quantun mechanical analogue of the second term of Eq. (2) would be of 
interest. In this context, a result first derived by Brow would prove to 
be useful. Brown 10 studied the K-shell asynptotic stopping power of an 
hydrogenic system (with 2 K-electrons) for a fast projectile, taking the 
maximum momentum transfer equal to 2mv- as if the electron was free. The 
asymptotic stopping power equation obtained by Brown can be expressed in a 
form similar to the Eq. (2). The first terms of both of these equations, 
since they involve the mean excitation energy, can be assumed essentially 
equivalent within the local plasma approximation. He obtained for the 
second term in Eq. (2) for an hydrogenic system, the quantity 1/n where, 

n = — m v 2 /Z 2 R where Z is the effective nuclear charge for s-sheil 
s 2 s s 

(s=K,L, ... ) and R is the Rydberg Constant. Walske, 11 on the other hand, 
took the upper limit for momentum transfer to be infinity, thus overestimat- 
ing the nuclear momentum recoiling and obtained instead 2/n g . Tn reality. 


however, due to the recoiling of the nucleus, the result should be expected 
to fall somewhere between 1/n and 2/n . This fact will be incorporated 
later (see Eq. (6)) as a parameter which we later estimate. At the 
present, however, for the sake of simplicity, combining Brown's result for 
the K-shell with Walske's result for the L-shell, i2 but retaining the 
consistency with the free election model, we vrite the analogous second term 
(known as shell correction) for an hydrogenic system with Z electrons as 


C = C K* total + C L 




) . 


(4) 


which can be rewritten for a real atom as 

C * h T“ ♦ (Z)', (5) 

“ A inv 2 

C 

where 4 (Z) = Z- f(Z)g (6) 

and <T> = - [zj*R + (Z-2) \ Z^R] . (7) 

In the above, a coefficient f{Z) has been introduced to distinguish a 
real atom from an hydrogenic one. The coefficient f(Z) is known to be 
less than unity for L-shells for low atomic nunber targets. The coefficient 
g is introduced to incorporate the effect due to the recoiling of the 
nucleus. 


6 



ji 


At this stage, it is appropriate to discuss various features associated 
with the low energy projectiles and the low atomic number targets. First, 
in the low energy region, the projectile's full charge z will not be oper- 
ational in the stopping process due to electron capture that is influenced 
mainly by the outer shell electrons of the medium. Second, Walske has 
pointed out that the coefficients f (2) are unreliable for the low atomic 
nunbers Z < 30 due to use of the hydrogenic wave functions. 

It is, evident from the above observations that some sort of crude 
estimate of the quantity C is in order. This is justified since the usual 
incorporation of these effects involves the fitting with the experimental 
data. The inclusion of the effect of projectile's effective charge should 
decrease the stopping number for all the elements; for Li the most and for 
Ne the least. In order to incorporate this effect and the other problem 
of the need for an accurate value of the coefficient f(Z) as stated above, 
it is reasonable as a first approximation to assume a semiempirical constant 
value of the quantity <j> (Z) equal to- half the total number of electrons in 
noble gas atoms. Such a division should over-estimate shell corrections for 
lithium and beryllium in decreasing fashion and under-estimate that for 
helium, neon, carbon, nitrogen, oxygen and fluorine also in a decreasing 
manner. Such a change in shell corrections is indeed \tfiat is needed as 
almost entirely compensating the effect of the effective charge of the 
projectile on the stopping power. Since, in this paper we are interested 
in the atoms with atomic number below 10 this assumption implies that 

♦ (Z) ■ 1 for Z < 2, (8) 

♦ (Z) - 5 for 3 < Z < 10. 




ji 

ii 
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Implicit in the above partition of <fi is the fact that the quantity C no 
longer represents the so-called shell corrections only, but presunably also 
some other effects including that due to the projectile's effective charge 
as well as due to the neglect of the higher order terms in Eq. (21, One can 
now write Eq. (5) as: 


and 


C = 1 <T> 1 fnv . 7 . o 

i f j-,‘i "* z « 2 - 

wmv‘ 


r 4 r, ‘ I for 3 < Z < 10, 


smv^ 


(9) 


where <T> by virial theorem is just the average kinetic energy of the 
electron and should be averaged over all the Z electrons in the atom. 

In order tc make a transition to an atomic system, we assune the above 
results and accordingly replace Eq. (2) with 


L 8 in Y 


L = An Y 




1 

7 


1 

Y 


for Z < 2, 


for 3 < Z < 10, 


( 10 ) 


The low and high energy L functions should now be combined to deter- 
mine the appropriate dependence of the stopping power on energy. To do 
this, we used Eqs. (3) and (10) for our desired results after replacing 


8 


by yojp in than, where nonconstant values of y were obtained from 
reference 2. Bonderup, 13 had combined Eqs, (2) and (3) and assuned a con- 
stant value of y equal to {f- Unlike Bonderup, we tried to preserve the 
continuity between the low energy stopping nunber function given by Eq. (3) 
and the high energy function given by Eq. (10). In this way, stopping 
number values for a system can be obtained given the velocity of the pro- 
jectile and the density p (?). 

For a diatomic molecule, the Gordon-Kim model gives the density as: 


p molecule ■ P a (r) + p fa (r -**>• 


(id 


where (r) is the atomic groundstate density. R ab is the distance 
between the two atoms, which is known to be 1.094A for N 2 and 1.207^ 
for the 0 2 molecule. Equation (11) was generalized for water vapor 
including its partial ionic bond nature and neglecting the overlap between 
the two H-atoms. The distance between the 0 and H nuclei was taken 
equal to 0.958A. 

The molecular stopping power for protons was obtained by averaging the 

•f 

stopping nunber over r for N 2 , 0 2 and water vapor molecules. Hartree- 

Fock wave functions were employed in these calculations. 

Table I lists the results of this paper, together with Andersen-Ziegler 
curve fitted results, 14 and two sets of experimental data for the 0 2 


molecules. Table II lists these values for the N 2 molecule. In table 
III, the results of this paper for water vapor are compared with the 
available experimental data for energy ranging from 40 keV to 500 keV. Good 


agreement, within 10 percent, is found with the experimental data. 

In order to discuss the departures from Bragg's rule, it would be 
relevant to cite a systematic study carried out in a series of experiments 
at Baylor University. 3 - 5 The study revealed that for low energy projectiles 
there may exist a deviation from Bragg's rule depending on the physical 
state, but most importantly, on the chemical structure of the compounds. 

The confusing status of the dependence on the chemical structure can best be 
described by citing these studies in chronological order. First in 1971, 
the Baylor group 3 sunmarized that the compounds with single bond and double 
bonds should obey Bragg's rule. The compounds containing triple-bond struc- 
ture were found to deviate from Bragg's rule by as much as 12. 8% (a p st- 
eles of energy between 0.3 and 2.0 MeV often were the projectiles). In 
particular, these authors indicated that the molecular hydrogen (single 
bonded molecule) should obey Bragg's rule. Later in 1972, the Baylor group 4 
critically looked again on their previous conclusions. They indicated that 
perhaps the hydrogen atomic stopping cross section may be considerably dif- 
ferent than one-half of the molecular stopping cross section and thus should 
cause considerable deviations. However, the Baylor group in 1974, 5 recog- 
nizing the difficulty of obtaining atomic cross sections exper imental ly, 
based their analysis on the existence of some modified, but unique atomic 
stopping cross sections. 

It is therefore imperative that in order to discuss the deviations from 
the Bragg's rule, one must have access to the atomic and molecular stopping 
cross sections. We calculated both the atomic and the molecular stopping 

cross sections as a function of projectile energy for the molecules O 2 , 

N 2 , and H 2 . These results together with the deviations from 8ragg‘s rule 


10 


are exhibited in tables IV, V, and VI. One sees that the deviations from 
Bragg's rule become small as incident energy increases — in agreement with 
observations made by many workers including those at Baylor University.’ It 
is to be noted that N 2 is a triple-bonded, 0 2 a double-bonded (from the 
bond energy point of view) and H 2 is a single-bonded molecule. The 
maximim deviations from Bragg's rule for energy 100 keV and above are 6 , 1 %, 
2.6%, and 10%, respectively, for these molecules. Thus, the deviation 
depends on the chemical structure. Wien the Gordon-Kim model is used, the 
overlap of electron density determines the deviation or molecular binding 
effects. For instance, for the hydrogen molecule, the distance between 
nucleons is very small, equal to 0. 74A. It is expected that the overlap of 
electron density is large, thus explaining the considerable deviation from 
the Bragg's rule. The stronger the bond energy, the shorter the distance 
will be. It is interesting to note that single-bonded, double-bonded and 
triple-bonded carbon molecules have internuclear distances equal to 2.94, 
2.52 and 2.24 in Bohr units, respectively. We may thus expect that the 
triple-bonded carbon will have more deviation from Bragg's rule than the 
single-bonded carbon. 
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TABLE I. Proton stopping cross section values (in units of 10" 15 eV cm 2 ) 
per atom. These were obtained by dividing by two, the stopping 
cross section per molecule of the oxygen molecule. 


E(keV) 

Theoretical 
Values of 
the Present 
Paper 

Curve 

Fitted 

Values 

Andersen 

and 

Zi egler 14 

40 

15.89 

14.6 

80 

17.48 

17.0 

100 

17.43 

17.0 

300 

11.84 

11.9 

500 

8.92 

8.8 

1037 

5.64 

- 

2591 

2.97 



EXPERIMENTAL RESULTS 

Reynolds 
et al . 1S 

Langley 16 

15.2 ±2.6 

- 

17.25±2.6 

- 

17.17*2.6 


11.99±1.7 

- 

8.84*1.7 

- 

- 

5.25 

- 

2.85 


n 
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TABLE III. Proton stopping cross section values (in units of 10" 15 eV cm 2 ) 
per molecule for water vapor. 


E(keV) 

40 

80 

100 

300 

500 

Theoretical 
Values of 
the Present 
Paper 

28.81 

27.8 

26.8 

17.1 

12.6 

Reynolds _et al , 14 

25.0±2.6 

27.6±2.6 

27.3±2.6 

17.9±1. 7 

13. 0±1 . 7 


O 
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TABLE IV. Deviations from Bragg's rule in the case of oxygen molecule. 
e(atomic), the atomic stopping cross section (in units of 
eV x 10“ l5 cm 2 /atom) was obtained from Eqs. (3) and (10). 

e (molecule) is in units of i eV x lO" 1 ? cm 2 /molecul e. 

2 


E(keV) 

40 

100 

200 

300 

500 

1037 

100000 

e( atomic) 

17.44 

17.48' 

14.65 

12.15 

9.1 

5.72 

0.1492 

e(molecule) 

15.89 

17.43 

14.36 

11.84 

8.92 

5.64 

0.1476 

percentage 

deviation 

8.9 

0.3 

2 

2.6 

2 

1.4 

1.1 


TABLE V. Deviations from Bragg's rule in the case of nitrogen molecule. 
e(atomic), the atomic stopping cross section (in units of 
eV x 10‘ 15 cm 2 /atom) was obtained from Eqs. (3) and (10). 
e (molecule) is in units of i. eV x 10~ 15 cm 2 /molecule. 


E(keV) 

40 

100 

200 

300 

500 

1037 

100000 

e( atomic) 

19.33 

18.57 

14.32 

11.56 

8.53 

5.30 

0.1340 

e (molecule) 

17.20 

17.79 . 

13.75 

10.85 

8.10 

5.20 

0.1319 

percentage 

deviation 

11 

4.2 

4.00 

6.1 

5.0 

1.9 

L3 





Deviations from Bragg's rule in the case of hydrogen molecule 
e(2xatomic), the atomic stopping cross section (in units of 
eV x 10" 15 cm 2 /atom) was obtained from Eqs. (3) and (10). 
e (molecule) is in units of eV x 10” 15 cm 2 /molecule. 


E(keV) 

100 

200 

300 

500 

800 

1037 

2591 

e(2 x atomic) 

12.7 

8.13 

6.1 

4.17 

2.89 

2.36 

1.11 

e(molecule) 

11.43 

7.53 

5.71 

3.93 

2.75 

2.24 

1.07 

percentage 

deviation 

10 

7.4 

6.4 

5.8 

4.8 

5.1 

3.6 


O 
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ABSTRACT 

A modified local plasma model based on Linhard-Winther, Bethe, Browi, 
and Walske is established. The Gordon-Kim's molecular charged density model 
is employed to obtain a formula to evaluate the stopping power of many 
useful molecular systems. The stopping power of and He gas was 
calculated for incident proton energy ranging from 100 keV to 2.5 MeV. The 
stopping power of 0 2 » N 2 and water vapor was also calculated for incident 
proton energy ranging from 40 keV to 2.5 MeV. Good agreement with experi- 
mental data was obtained. 

A discussion of molecular effects leading to departure from Bragg's 
rule are presented in this thesis. The equipartition rule and the effect of 
nuclear momentim recoiling in stopping power are also discussed in the 
appendix. The calculation! procedure presented in this thesis hopefully 
can easily be extended to include the most useful organic systems such as 
the molecules compound of carbon, nitrogen, hydrogen and oxyjen which are 
useful in radiation protection field. 
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INTRODUCTION 


The subject of energy loss of heavy ions such as protons, or a par- 
ticles passing through matter, has been studied for more than sixty years. 
Research in this area began with the study of a mechanism under which 
charged particles lose their energy mainly to the atomic electrons. These 
studies have contributed to the basic understanding of the interaction of 
charged particles with matter-atoms-molecules and, more recently, to materi- 
als. The energy loss parameters have found their use in various applica- 
tions. The list includes: radiation dose effects on solid state devices; 

shielding problem; space radiation research; design and calibration of 
spectrometers and dosimeters; proton doses in manned or unmanned space 
flights; energy transfer to living cells, and radiation effects in materi- 
als, etc. 

The Bethe [1] theory of energy loss of fast charged particles rests on 
the knowledge of the so-called mean excitation energy .of the medium. Once 
this parameter is known, the high energy stopping power of an atom can 
readily be calculated. The determination of this parameter, however, is 
very laborious, as is seen in the works of Dehmer, Inokuti, Saxon, and Baer 
[2], [3], [4], who calculated the mean excitation energy parameter for atoms 
of atomic number ranging from Z = 1 to Z = 38. The numerical evaluations 
involving the Hartree-Sl ater wave functions in these calculations are so in- 
volved that the estimating of the errors is difficult to ascertain. 

The Bethe theory, although developed for the atoms, has also been 
extended to obtain the stopping power of molecules under the Bragg's rule 
[5]. One essentially ignores the chemical binding of molecules under this 
rule. Recently, however, several experiments [6], [7], [8], have revealed 


that for low energy regions there may exist deviations from the Bragg's 
rule. Furthermore, there are some indications (see, for instance, a series 
of papers by Wilson and his co-workers) that the Bragg's rule may not be 
obeyed in the determination of the mean excitati. - energy parameter for 
molecules, although this departure does not have much effect on stopping 
power because of the dominance of the velocity of the projectile on the 
stopping power. 

It is evident from the above discussion that the traditional approach 
of obtaining the molecular stopping power from the atomic stopping power via 
the Bragg's rule should be abandoned, at least for the low energy projec- 
tile. The local plasma model which has been successful in predicting the 
mean excitation energy could serve as the appropriate candidate for an 
alternative approach. 

This report discusses the establishment of a modified local plasma 
model by employing the Gordon-Kim molecular [9] density model, Which 
provides a method of calculating molecular stopping power even at quite low 
proton energies. In spite of the fact that it is a somewhat average model, 
calculation is relatively simple, and the calculation of stopping powers of 
H 2 , He, 0 2 , N 2 and water vapor are in fair agreement with experimental 
data. Modifications to the local plasma model are mainly due to the com- 
plexity of real atomic and molecular situations. Besides the effects of 
shell corrections, and the screening of projectiles, other effects such as 
nuclear recoiling are also involved. Interesting discussions on deviation 
from the Bragg's rule and on the modified local plasma model and the conclu- 
sions deduced from experimental data are presented in the last section of 
this report, "Molecular Effect of Stopping Power." 

From the theoretical model established in this report, the departure of 
stopping power from Bragg's rule only occurs for low velocity projectile 


cases. Deviation from the Bragg's rule also is found to depend on the 
chemical structure of molecules. More overlap of electron clouds is found 
to cause more deviation from the Bragg's rule when the Gordon-Kim model is 
employed. The geometric structure gives the most important information on 
molecular effects. 

The basic stopping power theories assume the interaction between the 
projectile and the atomic electron, vrfvich is assuned initially to be at 
rest. While it is true that the overall energy loss must, indeed, take 
place in this manner, in some collisions the recoiling of the nucleus cannot 
be neglected. This is especially true for relatively low incident energies 
of the heavy ion. Although the recoil energy may be small due to the heavy 
mass of the target nucleus, the recoiling momentun may be very large. Thus, 
the conservation of momentun will lead to a different value of the momentum 
of electrons than has previously been assuned in these theories. 

This observation lends itself to solving a three body problem. Thus, 
in this report, a semi-classical three body model is established to calcu- 
late relevant quantities. Specifically, exact semi-classical three body 
calculations are made for a proton incident on a hydrogen atom with the 
electron in its first Bohr orbit. The model and the resulting conclusions 
especially in regards to the equipartition theorem for shell corrections are 
presented in the appendix. 

Historical Review 

The three prominent theories of penetration of charged particles in 
jj matter are: Bohr's semi-classical theory, Bethe's quantum theory, and 

Lindhard's local plasma theory based on the treatment of free electron gas. 

!j Since this paper touches upon all three theories, it is helpful to outline 
;! their main features and assumptions of interest. When appropriate, the 


detailed structure of these theories will also be given. In later sections. 


Semi-Classical Theory 

Bohr's semi -cl ass leal treatment of the slowing down charged particles, 
done as early as 1913 [10], [11], was the first to give the overall 
characteristic structure and the features of the penetration theory. The 
classical parameters appearing in this theory, surprisingly, could later be 
calculated or related to the quantum treatment. The main underlying 
assumptions and the characteristics of the theory are: 

1. The coulomb interaction between the incident-charged particle of 
velocity v and the atomic electrons is assumed to be responsible 
for energy loss. 

2. The momentum transfer is sufficiently small so that the projec- 
tile's path is a straight line. 

3. The atomic electron is ’assumed to be at rest. 

4. The minimum impact parameter is determined from the knowledge of 
the maximum momentum transfer of the projectile to a free electron. 
The maximum impact parameter is determined under the assumption 
that the interaction time must be larger than the orbital period of 
the atomic electron. 


Based on the above assumptions, Bohr obtained the stopping power formu- 
la as: 


£ 

dx 


4*NZ 


^e 4 

mv 2 



e 2 oi 


( 1 ) 


where, ze is the charge of the projectile and N is the number of atoms 
per unit volume, with Z electrons per atom, and u> is the characteristic 
atomic frequency. 

Equation (1), even though the derivation was based on classical consid- 


erations, embodies the main features of the quantun mechanical description 
given in the later section. As is known, in both theories, the properties 
of the incident ion such a? Its charge ze, and its velocity v occur in 
this equation, m being the mass of the electron. The properties of the 
medium are contained in the quantities N, Z, and w. 


Quantal Theory 

Bethe's non-relativistic calculation of stopping power was performed in 
1933 in the first Born approximation [1] and rests mainly on the following 
assumptions. 

1. The interaction responsible for energy loss is the coulomb inter- 
action between the incident ion and the atomic electron. 

2. The speed of the ion is much greater than that of the atomic elec- 
trons. 

3. The calculation of the maximum momentum transfer entails the colli- 
sion with the electron initially at rest. 

4. Within the plane wave Born approximation, the assumption that the 
electronic positions are correlated only over relatively small 
distances implies the use of the dipole oscillator strengths. 


Based on the above assumptions, B-'the arrived at the following energy- 
loss formula: 


d£ _ 4ir NZz 2 e 4 2mv 2 

dx mv 2 I 


( 2 ) 


where I, the mean excitation energy of the medium, is defined through the 
electric dipole strength f by: 


5 


■'VJcsTrr: 


jJ 


Z In I * 2 f In E_ 
n n n 


with E n being the ^enenergy of the electron. 

Notice the similarity between equations (1) and (2) obtained under 
classical and quantal considerations, respectively. Notice also the occur- 
ance of the factor 2mv 2 (which represents the maximum energy transfer of 
the assumption (3)) in the argument of log in equation (2). 


Theories Based on the Thomas-Fermi Model 


Bloch [13], in 1933, calculated the stopping power by using the Thomas- 


Fermi model for the many electrons of the atoms of the medium. Useful 


results were later obtained by Lindhard [14], [15], [16], and his co- 


workers. Lindhard showed that for a swift heavy particle of low charge the 


stopping power of a free electron gas is given by 


4 mz 2 e l+ 


. r 2mv 2 'i 
p in (— -) 

Yhu p 


where, p is the electron density, and w is the classical plasma 

o 4 7T e^ p 

frequency given by o»“ = p, and y is a parameter calculated by 

p m 


Lindhard to be equal to VT Lindhard made the bold assumption that the 


theory can be extended to the atoms if one whites the above equation as: 


4mz 2 e 1+ N 


/ In 


Y-hw (r) 


p(r) d 3 r 
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where, now p(r) is to be evaluated by employing the quant al wave functions 
of the electrons. This model is often called a local plasma model, in 
1 iterature. 

Further Discussion of the Bethe and Lindhard Theories . 

The simplicity in Bethe's theory results from several factors. The first 
factor allows all the atomic electrons to participate in the stopping pro- 
cess. Indeed, this is not valid for the inner shell electrons which are 
tightly bound and may not always participate in the stopping process. The 
second factor is the inclusion only of the dipole transitions in the theory, 
although other transitions (though less probable) may possibly take place. 
This, as indicated earlier, is tantamount to assuming that the electronic 
positions are correlated only over relatively small distances. These 
corrections to the Bethe theory have been extensively investigated. 

Lindhard 1 s local plasma theory, as mentioned earlier, makes the bold 
assunption through equation (3) for its application to atomic systems. This 
approach, surprisingly, works well for atoms as is shown recently in the 
evaluation of the mean excitation energy parameter [17]. The approach rests 
with the comparison of equation (3) of Lindhard with the Bethe formula 
equation (2): 

Z In I = / p(r) In (Yhu> p ) • d 3 r (4) 


The local plasma model is also relatively easy to extend to molecules for 


Ionic bonded gases, covalent bond gases and metals. Such progress has re- 
cently been made by various authors. (See, for instance, a series of papers 
by Wilson and his co-workers). 


STOPPING NUMBER FUNCTION 

The energy loss of fast charged particles caused by their inelastic 


collision with atoms, is given with good approximation by Bethe's formula, 
which is based on the quantum perturbation theory. However, the mean 
excitation energy in his formula is too complicated to evaluate theoretic- 
ally for many practical problems especially for molecules. It appears the 
local plasma model affords a simple method to calculate the mean excitation 
energy of elements as well as compounds. There still remain three problems: 
1) people have more interest in stopping power than in mean excitation 
energy; 2) how to extend this method to evaluate molecular stopping power; 

3) how to extend this method to slow charged particle cases as well as to 
fast charged particles cases. Now let us recall the basic formula of local 
plasma model for fast charged particles. 


dE 

dx 


4 wz 2 e 4 N 

mv 2 


\p(r) d 3 r 
Yt*»p(r>y 


(5) 


2 — Ck 2 q ( \ 

where ^ p (r) = — v ' is the plasma frequency. We can see in the above 
formula that once the parameter y and the electron density p(r) are 
known, the stopping power can be determined immediately. We will discuss 
the problem of parameter y in the section "Local Plasma Model for High 
Projectile's Energy Cases." 

The atomic electron density is just the square of the wave function of 
the ground state of the atom. It is not easy to find the molecular wave 
function. Fortunately, we have a very simple molecular electron density 
model, namely the Gordon-Kim model [9]. It is a rough model, but as a first 


order approximation this model gives the main feature,:* of molecular effects. 
Detailed discussion will be presented in the next chapter. To extend this 
model to low energy regions is of considerable interest from both a practi- 
cal and theoretical point of view. Unfortunately, many effects arise at low 
energies including projectile charge screening, nuclear momentun recoiling, 
forbidden transition, etc. On one hand we have to consider the above 
complexity in realistic atomic and molecular world. On the other hand we 
still need to keep the simplicity of the local plasma model, otherwise no 
j results can be obtained in practice. This section is devoted to establish a 
modified local plasma model. We will concentrate on finding a stopping 

K 

, nunber function L. It should be valid for low energy as well as high energy 
cases, and should still retain the simplicity of the local plasma model. It 
l will also approach the realistic cases as closely as possible. However, it 
’ is only an average model and some estimating is involved. 


Local Plasma Model for High Projectiles Energy Case 
When the projectiles move rather fast the stopping power can be deter 
mined by Bethe‘s formula with good accuracy. 



where I, the mean excitation energy has been of considerable study. This 
l is because once the parameter is known, the energy loss can be obtained 

ji 

immediately. However, the evaluation of this parameter is considerably 

\\ 

5 

! laborious. In principle, the parameter can be calculated exactly only for 


the H-atom for which accurate wave functions are known. More recently, 
Dehmer, Inokuti and Saxon [2] used the tabulation of Hartree-Slater poten- 
tial given by Herman and Ski liman. They solved the Schroedinger equation to 
obtain the radial matrix element R (nl, n'l') where (n, 1), (n‘, 1') are 
initial and final states of atom. Knowing these values they calculated the 
mean excitation energy parameter I for atoms ranging from Z=1 to Z=38. 
However, there are difficulties with the above approach. In practice, the 
extreme complexity of nunerical calculations renders it impossible to extend 
the approach to evaluate the mean excitation energy of molecules. 

An alternative approach is the possibility of using the local plasna 
model, formula (5). There are, however, two central quantities which 
should be known. These are function y and p(r). The evaluation of p(r) 
rests on the determination of the wave functions of electron' only in ground 
states, such Hartree wave function for various atoms have recently been 
available in the work of Clementi and Roetti. The problem is then the 
determination of parameter y. 

Calculation of Parameter y 

In his original theory of stopping power of electron gas Lindhard gave 
a quantitative discussion for the values of parameter y. The first term in 

<pmy2 

high energy expansion of stopping nunber L is given by In C(x)) 

where C(x) =1. . Linhard discussed the function C(x) as a function of 
electron density. He surmised that the function C(x) depends on the den- 
sity slightly only. It should be a little bit less than one for moderate 
density and should approach unity for both extrG’nely low and extremely high 


densities. Furthermore, Lindhard suggested that the values of y can be 
taken to be unity for light atoms without large error. Later, Lindhard and 
Scharff [14] suggested on the basis of a simple model that y = vr can be 
chosen for heavier elements. 

Chu and Power using constant value of y =VT obtained the parameter 
I for various atoms. [17] Their calculated values however were found to 
exceed the values calculated by Delmer et al. (based on oscillator strength 
method) by 20% to 30%. The most satisfying thing, however, was the similar 
trend in the variation of I values in both cases as a function of atomic 
nunbers. This points towards a greater confidence in the local plasma 
theory. Some authors encouraged by this, and also not satisfied with Chu 
and Power's results treated the parameter y as an empirical parameter to 
fit the data. Unsatisfied with this type of empirical treatment, we adopt 
the following approach. Bohm and Pine [18], [19] have treated the problem 
of collective long range interaction in a quantun electron gas. They 
introduced normal coordinates of collective motion of electron beyond some 
screening radius r c and the individual particle motion was considered to 
be more important for radius smaller than r . 

Bohn and Pine showed that the average plasma frequency <w> (which is 
the average over the frequency of the collective oscillations, say w) is a 
linear function of the classical plasma frequency o) p and is given by 

<o>> = [1 + 1 £ (1 + -2 8 2 )] a) (7) 

2 X_ 10 p 

s 

where X. = r /a is a dimensionless parameter and is the average distance 
s s o 


between electrons, and where a 0 is the Bohr radius. 


The parameter 0 can be determined by minimizing the electron long 
range correlation energy {the long range part of correlation energy is 
obtained by subtracting short range exchange energy from the cohesive 
energy) 


E 


lr 

corr 


0.866 e 3 

,3/2 

s 


0.458 0 2 + 0.019 8 4 


( 8 ) 


The minimization of the above equation leads to the following equation for 3 


0.076 0? 


+ 2 - 598 3 



- 0.916 = 0 


(9) 


At this stage let us recall that Bohr's semi-classical theory and Bethe's 
quantum theory have the similar first term in the high velocity expansion 
of 

L = In (2SlL) 

One must note that E is some average energy which in the case of the Bethe 
theory is the mean excitation energy I. Thus, some sort of average quantity 
must occur in the argument of log term, therefore, recalling equation (4) 
of Lindhard' s theory 


It was show by Bohm and Pines that yfw may be replaced by <w> where 

g2 

<w> is given by equation (7). Thus y = 1 + 3/2 (1 + 3/10 0 2 ). Once we 

a 

know A then 0 can be determined by equation (9), hence y can be 
•# 

determined by equation (7)* Now let us use this average model to evaluate 

y as function of Z for various atoms. We took average distance between 

3 4 3 

electrons in an atom as r $ = ro/Z where ro is atom's radius. Then 
X s = r s /a o w * 1ere a o is So ^ r rad i us . Table 1 shows r , $ , x $ as 
functions of Z for some selected values of Z ranging from Z=2 to 


Table 1. r g , 8 as 


I 


z 

r s 

A s 

2 

1 

2.42 

3 

1.55 

3.27 

4 

1.12 

2.15 

5 

0.98 

1.75 

6 

0.91 

1.53 

7 

0.9 

1.43 

10 

1,17 

1.65 

11 

1.9 

2.60 

12 

1.6 

2.13 

13 

1.43 

1.85 

14 

1.32 

1.67 

15 

1.28 

1.58 

16 

1.27 

1.54 

18 

1.43 

1.66 

36 

1.59 

1.47 

54 

1.75 

1.41 


| Table 2 shows y as a function of Z. From table I 
f a constant equal to 1.19. This is reasonable, si? 


chosen values of y ranging from y = 1.1 to y = 


8 


0.54 

0.62 

0.51 

0.46 

0.43 

0.42 

0.45 

0.55 

0.50 

0.47 

0.45 

0.44 

0.43 

0.45 

0.42 

0.41 











Lindhard's description of y it is a function that slightly depends on 
density or X , and is not a constant. We obtained y * 1.19 partly 
because we took an average model. Table 3 lists y as a function of X , 

*5 

(ranging from 10- 18 to 10 7 ). Since p(r) * -i-- from the table we can see 

x s 

that y is very slightly dependent on density; a little bit greater than 
one for moderate density and approaches unity for both extreme cases - very 
low and very high density. Compared with Lindhard's qualitative description 

of C(X) = in his theory, as mentioned earlier in this section, we 
reach agreement with Lindhard for these two different cases. In the exact 
calculations of mean excitation energy using local plasma model, we took y 
as fvCv,Lion of distance r (distance from electron to nuclei), instead of a 
constant. 

Recent development confirmed the concept of y obtained here in 
improving the work of Chu and Powers. Furthermore, the above consideration 
will be used in the following section for low energy stopping power. 

Low Energy Stopping Number Function 

Lindhard-Winther expanded the stopping number function of free electron 
gas in high energy case and in low energy case as follows: 

L 2 * In f 2mv ) - <T> f m (high energy case) (10) 

V / 2 
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Table 2. y as functions of 





















Li 


(— ) 7 C,(X) (iffiii) 7 
3 h- p 


(low energy case) 


where Cj(X) 



[In ( 


2 

1 +7 X 2 
X 2 


) - 


X* 

1 - T 


] 


(ii) 


( 12 ) 


X 2 * 


TrhVp 


density V, 


* w p = "m ' 6 ’ * * 1S P^ asma frequency, p is the electron 

is the Fermi vel j ity. 


Bonderup [20] directly combined low energy L function and high 
energy L function L 2 and used VT instead C(X) in the first term in 
I 2 • (In [15] Lindhard and Winther assumed C(X) = 1 see last few lines of 
page 10). He then performed calculations for the stopping power of some 
elements using the above L function and local plasma model. Good agree- 
ment was found with experimental data for proton energy over 500 kev. Un- 
fortunately, for low energy regions, his simple approach is not valid. 

We establish our modified L function based on the following 
principl es: 

1. smoothly join Linhard-Winther' s high energy L function L 2 and 
low energy function L^; 

2. involve Pine correction using yha» p instead of h<Dp in all the 

terms in L function; K M 

3. apply correction on the second term of L 2 . 
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The correction on the second term warrants some discussion. 

Many complicated physical effects are involved in the low energy region 
We will discuss them in detail to approach the realistic atomic and molecu- 
lar world. Bethe, Walske, and Brown [21]* [22], [23], developed quantun 
mechanical theory of stopping power of innershells of atoms. Calculations 
are done under plane wave Born approximation using hydrogenic wave func- 
tions. Brown [21] and Walske [22] had calculated the stopping power of K 
shell electron. Walske [23] also calculated the stopping power of L shell 
electrons. Khandelwal and Merzbacher [24] calculated the stopping power of 
M-shell electrons. Khandelwal more recently evaluated K and L shell 
corrections [25]. 

They defined a stopping nunber function B as follows. 


dE 4ire 4 z 2 u 
- — * N z 

dx mv 2 i=K,L,M 


(13) 


where B, = / W| " ax wdw A 9 * |Fw, (q) | 2 1*1 (14) 

w m1n %ii n ** 

where | Fw- ( q) | 2 is the form factor, q is the change in incident parti- 

1 

cles momentun divided by (2m Z? eff . Ry) . For high energy projectiles, 
they also give an asymptotic stopping nunber formula of the form 
B s <e s , n s ) = S $ (0 s ) In n $ + T s (e $ ) - C $ (e $ , n $ ) where s = K,L,M and 
where C $ are the so called shell correction terms. Equation (13) can also 
be written in terms of 


where I is the mean excitation energy, e $ is dimensionless screening 
parameter denoting the observed ionization potential in units of ideal 
ionization potential Z 2 sRy/s 2 . n is also a dimensionless quantity and is 
given by L mv 2 divided by Ry. 

The calculations of the shell corrections are extremely complex. 

Bethe, Brown, and Walske expanded them in power series of . The first 

n s 

term of shell correction given by Brown [21] is C u total ~ for K shell 

K n k 

2 

electron meanwhile Walske' s [23] results are C. total ~ — The difference 

k n L 

between Brow and Walske is mainly due to the fact that they took different 
upper limits to estimate maximum momentun and energy transfer. As is well 
known when a heavy particle collides with a free electron, the maximum 
momentun transferred to the electron is 2mv and maximum energy transfer is 
2m v 2 these values are the ones taken in by Brown in doing his calculations. 
But if one considers the nuclear recoiling and the binding effects, then the 
upper limit of momentum transfer is no longer 2mv. Walske took the upper 
limit of both momentun transfer and energy transfer as infinity. Indeed, 
both Brown and Walske simplified the problem in this manner. Later, we'll 
use an exact three body semi-classical model to estimate the upper limit of 
momentun transfer. Right now, at this stage, we'll use Brown's assumption 
for consistency with the free electron gas of plasma model, but keep in mind 
that there are some errors due to Brown's assumption. The accurate result 
should be expected to be between Brown's and Walske’ s results. Now we try 
to establish some relationship between plasma model and quantum mechanical 
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calculation based on Bethe, Brown, and Walske's theory. Following Brown 


'k total ’ 


i _ 1 <n> 


n k <n> n k 


where we define <n> »— mv 2 /<T>, <T> Is the average kinetic energy of 

2 

electron of atoms then 


1 


<n> , 
n,. 


"w 2 


<T> 


keff 
1_ mv 2 
2 


by virial theorem <T> = (<E>j we have 


<n> 

n^ 


Z keff 

“F E5T“ 


For Z = 1 we get 


hence, we have 


^klz=l 2 ^ total 


C| Jz=i 

Z=1 2 <n> mv 2 2 


a i — L. = = i (Li nhard-Wi nther' s second term of L 2 function) 


For Z = 2 we get 


C kU=2 " C k total 


so that 
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^z*2 .1 JL_ 
Z, z=2 2 < n > 


1 

2 


(Linhard - Wi other's second term of L 2 


function) 


Now, let us consider atoms which contain L shell electrons. 

Walske also showed that the total shell correction of L electron 
C L total ~ njj ■ which ’ s still due to the value infinity for the upper limit 

of momentun transfer and energy transfer. If we keep Brown's model, i.e. 

assune that the incident particle collides with free electron, we have 

1 2 

f 1 , 7 mv 2 

L total ~ n, where n L = ; — • Now consider L shell closed. 

L Z Leff 2R * 

The total K and L shell corrections is 


: K total + C L total = “ + 


= l_ 


[1 + ] , _L [14. Z Leff 2 i <n> 

Z Keff 2 <n> " 


’Keff 2 


<n> 

n, 


= _L_ [l + Zl - eff2 i 2^ mv 2 /<T> 

<n> Z o 1 

KSff2 r^Wky 


= n - 4 - Z leff 2 1 Z Keff 2 Ry 

<n> z „ 

Keff 2 


|<E>| 


for h)drogenic atom 
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|<E>| * 


22 


Keff 2 


Ry + 4 Z Leff2 Ry • 8 


2Z Keff zRy * (1 + — ff2 - ) i 

W z 


hence, K total + L total _ _1 1 

l 2 <n> 


Now consider L shell open: 


C K total + C L“ - *L -L[l ♦iiHA] 

\ n L 8 n K 8 \ 


= J_ [1 + (Z-2) Z Leff 2 i Z Keff 2 Ry 


<n> 


8 Z 


Keff 2 


|<E>| 


where. 


. z Keff 2 «y _ <jj> 

l«>l ’ n„ 


N°« |<E>| - J [2Z Keff2 Ry + i. Z Leff2 Ry (Z-2)] 


hence. 


* |Wf zia + Zf mfi 

8 Z Keff 2 
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c k total + C L 


. J_ [1 +1111 Z| - eff2 ] 1/ M + (Z-2) Z Leff 2 ) 
<n> 8 Zfceff 2 2 8 Z Keff2 


„ C K + C L . 1 1 

or * — 

Z 2 <n> 

In both cases either closed shell or open shell we obtained the sane results 


C 1 1 1 <T> 

7 “ 7 <n> ”2" T 
L 2 mv 2 


= i (Li ndhard-Wi nther' s second term 
mv 2 of L, function) 


(15) 


It appears that from quantum perturbation theory we obtained the 
similar result as of plasma theory except a factor of ~ . Actually if we 
use Walske's result exactly the sane results are obtained from both quite 
different approaches. What is implied in this surprising similarity? We 
believe this similarity is the real background of local plasma model. In 
plasma model average kinetic energy <T> is related to the plasma density. 
On another side in Bethe-Brown-Walske' s theory as we have shown the average 
kinetic energy <T> is related to atomic wave function by using virial 
theorem or averge kinetic energy is related to atomic electron density. 
Since these two approaches 'give the same stopping power expressions, we may 
say that plasma model can be localized by the equivalence of these two 
approaches. But for low energy projectile, especially for light elements 
which contain L shell electrons, this formula is too simple to describe 
the stopping power due to various physical effects. First, one should take 


into account the transitions that are forbidden by Pauli's principle (for 

the importance of this in the case of asymptotic stopping niwbers see 

Khandelwal [26]). There is a modification factor on — or i- or — 

n K \ n M 

which depends on Z and also depends on the shell. The following tables 4 
and 5 list the coefficients of n L for different (where e L is the 
energy difference between ground state and lowest occupied state in units of 
\ e ff 2 ^y^ f° r L shell n = 2) according to Walske for 0 L = 1 the 
coefficient of n^" 1 is 2. 


Table 4. Coefficent of n^" 1 


S L 

0.35 

0.45 

0.55 

0.65 

1 

Coefficients 
of n L “ l 

1.5032 

1.0756 

1.9890 

2.0000 

2 


Table 5. Coefficent of 


0 L 

0.7 

0.75 

0.8 

0.85 

0.9 

Coefficients 

of V 1 

2.0662 

2.0999 

2.1196 

2.1290 

2.1309 
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Second, at low energy we also need to consider high order expansion 
terms which are not completely known. Walske's expansion for the L shell 
for the n^’ 4 forms involves fitting. 

Third, Brown and Walske's calculations are based on hydrogenic wave 
functions, as Walske has pointed out, for ZOO the results are not 
accurate. 

Fourth, as mentioned. Brown's model of upper limit of momentim transfer 
was used and this involves an uncertain factor between 1 and 2. This is 
because due to the nuclear recoiling, the real value should be between that 
of Brown and Walske, and should depend on projectiles energy. 

Furthermore, there exists a screening effect on low energy projectiles 
charge. It is easy for the element like Li to loose its valence electron, 
when positive heavy ion projectile moves rather slow. There is some chance 
that the projectile can capture the eletron. This will cause screening 
effect on projectile which in turn will decrease the stopping power. All 
these complicated effects should be taken into account. 

In formul a (15) 

C K + C L _ 1 r 1 Z, 

— 1 7 t<n> T 

s Z {j Li nhard-Winther' s Second term of L 2 function times Z} 


1 C 

the factor j is due to the definition of j as another Z in the 
bracket is mainly due to each electron's contribution to total shell 
correction, “—or is full shell's shell correction, it looks like 


or 


these two Z factor cancel each other but taking into account all the above 
effects the situation is not so simple. 

As we know the quantity Z in the bracket of formul a (15) is mainly 
due to each electrons contribution to total shell correction but at least 
the screening effect has a negative influence. As is well known, Li's 
first ionization potential is about 5 ev. It means only 5 ev additional 
energy c.an cause Li to loose its valence electron. Meanwhile, neon's first 
ionization potential is about 20 ev. The ionization potential rises linear- 
ly in Z between lithiun and neon. It appears for Li screening effect 
is much more important than for Neon. At least for 2 < Z < 10 elements of 
small Z have more influence on reducing the stopping power through screen- 

j,- ing than large Z. 

t, 

‘ Formally, we absorb this effect into the bracket. Now, in the bracket 

u 

we have two contradictory effects. The shell correction proportional to Z, 
the screening effect is in the inverse direction. In addition, there are 
many other unclear effects which we have mentioned above. Roughly now, we 
have the following assunption, suppose these two contradictory effects 
roughly cancel each other also with other effects. We may use a half full 
shell nunber as an average instead of Z in the bracket. 

Thus, we obtain the expression for L 2 : 


in - i 52. for 2 < 2 

hcu Z mv 2 
P 


( 16 ) 


u • in - £ £1 


for 2 < Z < 10 


where <T> = i 
mv 2 2 


(Lindhard-Winther' s second term in L 2 function) 


As mentioned earlier the equivalance of the shell correction term from 
Li ndhard-Winther' s theory and Bethe-Walske' s theory actually gives somewhat 
the explanation of the local plasma model. Furthermore, for real atoms we 
also give the correction on Lindhard's theory. Also we can tell how far the 
local plasma model can be applied to real situations. Lindhard and Winther 
also used a parameter y in both the first term and terms of L 2 to 
simplify the calculations we will show in the note at the end of this sec~ 
tion in detail. 

Now (16) becomes 

L 2 * In Y - i - i for Z < 2 

Z 10X Y 


1-2 




1 

Y 


for 2 < Z < 10 


where X 2 = f_ 
irhV,. 


Vp is the Fermi velocity and Y = 


2mv 2 

hw 


Now we apply the Pine correction on L, i.e. use yhco p instead of hu» p 

2mv 2 ?mv 2 

in L 2 as well as in L l or using Y = — — instead of Y = in L 


Yhw 


hw. 


function. 

\ 

Finally, we smoothly join L l and Lj thus obtaining L function by 
the following manner. Extensive nuneri cal evaluations using a computer 


program of the function Li and L 2 for various values of the variable Y 

revealed that in most useful cases there were found two roots Yi and Y2 
of the equation Li = L 2 where Yi < Y 2 . Furthermore, the slope M of 

the function L2 for values of Y greater than Y2 was always small. On 
the other hand, the slope of L 2 for Y values less than Y A was very 
steep. These observations, including the behavior of these functions (see 
figure 1 ), led us to the following recipe to preserve the continuity consis- 
tent with the physics of the situation 


L * L2 when MY < L2 

L 2 = Max [Li, L 2 ] when MY > L 2 


actually we have 


Li V < Yi 

L 2 Y > Yi 


vrtiich was found to be convenient for the conputer programming. 
The functions L A and L 2 are defined as follows 


3 3 

X 2 4 2 

Li = (L.)' t Ci (X) 
3 


1 


1 + ix 2 
3 


where Ci (X) • = — « — (In ( 2 — ) - 


X 2 « 
2 ( 1 - — — ) 2 

3 


X 2 


1 - ± X 2 
3 

1 + £ X 2 
3 


) 
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* 


l 2 * in y - 1 £1 i 
Z 10X X 

where N is the half full shell nunber 


N = 1 Z < 2 

J 5 for 2 < Z < 10 

NOTE: Lindliard and Winther used a parameter Y in both first and second 

terms of the L 2 function, i.e. 


instead of 


L 2 * In y * 


3 1 " 5 1 
5x Y 


U = In 


2mv 2 <T> 

hui 1 2 

p — mv^ 
2 


t/ 


where Y = 2myZ and X 2 = e2 


hw 


*hV, 


where V p is the Fermi velocity V P = h (3ir 2 n) 3 and w 2 

h m P 

n is the density of plasma. 


4Tre‘-n 

m 


i ; * 


* * 

i t • 
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Applying this formula in the atomic or molecular scale, i.e. in the 
local plasma model n is the electronic wave function square of ground 
state. Thus Y and X are the function of density since we know the wave 

function, L 2 can be easily obtained. Now we will show that the second 
term of l 2 

r-^— — is equal to - i 

- mv 2 5 X Y 

2 

the following formulae are used 


X 2 * - — definition of X 

*hVp 


O). 


2 _ 4rre 2 n 


m 


e_ = i m V P 2 = — (3ir 2 n)^ 
F 2 F 2m 


<T> = _2 m V P 2 
10 *■ 


From (17.3) we have Jn~= Vp •JvJT 

™ tt h T 3h ’ 


from (17.2) we have u> 2 = 

P 


4*e 2 


m 


(17.1) 


(17.2) 


(17.3) 


(17.4) 
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h “ p m y 

'r-s Vs ’'Vv --V5T 

where A * 2 

me 4/— 
y 3nh 

v . w 

h “p 

VvT- 2m ^" <» )* 

A Vp2 ^ F A Vp 

from (17.1) we have = A .J ®i 

X TT h 


thus 


Y -- V3~ ^ 

X V_ 


since from (17.4) we know <T> = -1 m , finally we have 


U * In Y - 


,1.5 


5x Y 


rm 


t 
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CALCULATION OF MOLECULAR STOPPING POWER 
In this section, we perform calculations on the molecular stopping 
power of H 2 , He, N 2 , O 2 and water vapor for wide proton projectile energy 
regions (for H^,, He, 100 Kev - 2.5 Mev, 40 Kev - 2.5 Mev else). Comparisons 
of the experimental data are indicated in tables 6 to 11, (see reference 
[27] and [28]). 

The basic formula for the stopping power of local plasma is as follows: 

. £ . j p(r) d 3 r 

dx mVp 2 Yhw^r) 

To extend the formula to low energy regions, it is only necessary to replace 
2mv 2 

In ( J- — ) by L function, which we have established in the section en- 

Yho> p (r) 

titled "Stopping Number Function." 

Now we are interested in molecular stopping power. The charge density 
in the integral should be the electron density of molecules. As we mention- 
ed in the section titled "Stopping Number Function," the Gordon-Kim model 
is a rough model of molecular electron density. They assuned that no 
arrangement or distortion of the separate atomic density takes place when 
the atoms are brought together. The total electron density of two 
interacting atoms is, therefore, taken simply as the sun of the two atomic 
densities. 

P (r) = p(r) + p (r - R 12 ) (18) 

molecule 


where Rio is the internuclear distance. 
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Let us use H 2 as an example to explain how to apply formula (18). 

1 — p 

The ground state wave function of the hydrogen atom is — e . The elec 

ir 

1 -2r 

tron density of hydrogen atom is — e . According to (18) 

IT 


, 1 „-2r 


H2 


_ e 


I 2 * 

-(^/r 2 +fti 2 - 2 rRi 2 cos^) 


V 


For other atoms the atomic wave functions are obtained by Clementi, Roetti 
[29] (published in atomic data and nuclear data). They used Roothan- 
Hartree-Fock method to calculate basic function and their coefficients for 
ground and certain excited states of neutral and ionic atoms for Z< 54. For 
molecules which contain more than two atoms the model can be extended as 
fol lows 

p molecule = P A^ + P B ^BA^ + p C^ r " *AC ^ 


Here we present the calculation of water vapor which contains three atoms. 


The angle between two OH bounds is about 105*. It should be calculated by 
formula (19), but for simplicity in this section we neglected the H-H over- 
1 ap simply as fol lows 

P H 2 0 = P o (?) + 2 P H (? - 0R) 

Also we need to consider partial ionic bond effect. Many compounds exhibit 

partial ionic bond rather than pure covalent bond. Pauling defined partial 

ionic fraction [30] as the measured dipole moment divided by ideal dipole 

moment. Consider, for example, the compound, HC1 . The distance between two 

nuclei is 1.27EA*, the partial ionic fraction is p = ^-«x 100%, y * = e x Ro 

* 4.80 x 1.275 = 6.12 D, but the actual measured dipole moment y* = 1.03 D. 

1 03 

Hence, p = x 100% » 17%. HjO is also partial ionic bonded compound. 
The actual measured dipoment is y = 1.940 the ideal dipoment y° = e Rq = 
4.8 x 0.958 D = 4.5980 the partial ionic fraction p = — = 0.42. The 

M * 

mole-cule's electron density of partial ionic bonded compounds can be 
expressed as follows 

p molecule = P A^) + P B ^ 

+ — + _ 
with p~= (1- p) p (r) + p p (r) 

The case p = 1 (pure ionic bond) physically corresponds to moving one 
electron from one atom to another atom. Now for one atom we use wave 


function of positive ion, for another atom we use wave function of negative 
ion. For the case p=0, that is the pure covalent bond, we still use the 
neutral atom's wave function. Hj, 0 2 , N 2 » are all pure covalent bond 
molecules. 

Water vapor is a partial ionic bonded compound. The distances between 
nuclei are listed in table 6. 

Table 6. 



H-H 

H-0 

— 

N-N 

0-0 

R(A) 

0.74 

0.958 

1.094 

1.207 


Table 7 lists the results of this paper together with Andersen and Ziegler 
[31], curve fitted results and three sets, of experimental data for H 2 
molecules. Good agreement within 20 percent is found with experimental data 
from proton energy 100 Kev. - 2.5 Mev. Calculations from the equations 
established by Bonderup [20] were also undertaken by extending them to 
molecular system. Table 7 lists these values for H 2 molecules in the last 
column. These differ from our results in the low energy regions. Table 8 
lists the same physical quantities for helium gas. The same trend is 
observed as in the case of the H 2 molecule. 

Table 9 lists the results of this model together with Andersen-Zeigler 
and three sets of experimental data from proton energy 40 Kev - 2.5 Mev for 



N 2 molecule. Table 10 lists the same physical quantities for 0 2 . Table 11 
lists the results of this model together with Reynolds et al. experimental 
data for water vapor all these tables show theoretical results of this 
model are in good agreement with experimental data. 
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♦Andersen and Ziegler [31] note that the data in reference [32] appear to be 
constantly high by amount ranging from 5 percent to 10 percent. 













m stopping cross section values (in units of 10” 15 eV cm 2 ) 
ire exactly A. stopping cross section per molecule 

5 eV-cm 2 ) of^N 2 gas 


E(keV) 

Theoretical 
Value of 
the Present 
Paper 

Curve 

Fitted 

Values 

Andersen 

and 

Ziegler [31] 

Best Avail. 
EXPERIMENT/ 

able 

M RESULTS 

t — - ' 

Reynolds 
et al. 
[32] 

Philips 

[36] 

Langley 

[33] 

40 

17.20 

16 

17.1 ±2.6 

14.1 


50 

17.81 

16.9 

17.8 ±2.6 

14.8 


60 

18.24 

17.3 

18,2 ±2.6 

15 


70 

18.48 

17.8 

18,5 ±2.6 

14.9 


80 

18.41 

17.9 

18.5 ±2.6 



90 



18.25±2.6 



100 

17.79 

17.7 

17.9 ±2.6 



200 

13.26 

14.1 

14.7 ±2.6 



300 

10.85 

11.2 

' 11.2 ±1.7 



400 

9.24 

9.3 

9.34±1. 7 



500 

8.10 

8.1 

8 . 08±1 . 7 



600 

7.25 

7.2 

7.21±1. 7 



700 


6.7 




800 


6.0 




900 


5.65 




1037 

5.20 




4.78 

2074 






2591 

2.71 

2.72 












Table 10. Proton stopping cross section per atom (in units of 10" 15 eV cm 2 ) 
or more exactly! stopping cross section per molecule 

(IQ-* 5 eV“cm 2 ) o? 0 2 gas 


E(keV) 

Theoretical 
Value of 
the Present 
Paper 

Curve 

Fitted 

Values 

Andersen 

and 

Ziegler [31] 

Best Av 
EXPERIME 

iil able 
:ntal RESUL1 

rs 

Reynolds 
et al 
[32] 


Langley 

[33] 

40 

15.89 

14.6 

15.2 ±2.6 

12.5 


50 

16.52 

15.5 

16.4 ±2.6 



60 

16.99 

16.2 

16.9 ±2.6 

14.2 


70 

17.29 

16.7 

17. 1 5±2 . 6 

13.8 


80 

17.48 

17 

17.25±2. 6 

13.8 


90 


17.1 

17.25±2.6 



100 

17.43 

17.0 

17.17±2.6 



200 

14.36 

14.6 

14.7 ±2.6 



300 

11.84 

11.9 

U.99±1.7 



40C 

10.14 

10 

9.76±1. ? 



500 

8.92 

8.8 

8.84+1.7 



600 

7.99 

7.9 

7.9U1.7 



700 


7.0 




800 


6.5 




900 


6 




1037 

5.64 




5.25 

2074 






2591 

2.97 




2.85 






« 

Table. 11. Proton stopping cross section per molecule (10" ls eV"cm 2 ) 
of HjO vapor 


E (keV) 

Theoretical 
Values of 
Present Report 

Reynolds 
et al. [32] 

40 

28.81 

25.0 ± 2.6 

50 

28.81 

26.1 ± 2.6 

60 

28.59 

26.9 ± 2.6 

70 

28.22 

27.5 ± 2.6 

80 

27.77 

27.6 ± 2.6 

90 

28.28 

27.5 ± 2.6 

100 

26.77 

27.3 ± 2.6 

200 

21.04 

22.0 ± 1.7 

300 

17.06 

17.9 ± 1.7 

400 

14.43 

15.0 ± 1.7 

500 

12.59 

13.0 ± 1.7 

600 

11.20 


700 

10.13 


800 

9.28 


900 

8.56 


1000 

7.97 
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MOLECULAR EFFECT OF STOPPING POWER 
Deviations from Bragg's rule due to chemical structure have been 
recently systematical ly studied by several experimenters. Soma interesting 
conclusions were summarized from these experiments. Now, we have already 
obtained stopping cross section of molecules by using method established in 
the sections titled "Stopping Nunber Function" and “Calculation of Molecular 
Stopping Power.” We can also easily obtain the atomic stopping cross sec- 
tion by calculating the difference of these two quantities. In this way we 
can obtain the deviation from Bragg's rule theoretical ly. A discussion, 
based on this evaluation and on conclusion from the experimental side, will 
be presented in this section. 

Discussion on deviation from Bragg's rule due to chenical structure 
of stopping power of molecules is of great interest. In 1905, Bragg and 
Kleemman first proposed the Bragg's rule [5 ] . It states that the stopping 
power (or stopping cross section) of a molecular substance is the additive 
sun of the atomic stopping powers multiplied by the nunber of times each 
atom occurs in the molecule. Bragg's rule has been shown by Thompson [38], 
for very high velocity proton to be valid within about 1%. Wilson and his 
co-workers [39], [40], [41], [42] systematically studied the molecular 
effect on mean excitation energy. Considerable deviations from Bragg's rule 
on mean excitation energy were found. This is not contradictory to Thompson 
because due to Bethe's formula 

dE 4rre 4 Z ^ ^ n 2mv 2 

dx mv 2 I 
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even though there may be considerable difference in mean excitation energy 
I, for very high velocity case the percentage deviation on stopping power 
can still be very small. 


For low energy projectile the situation is more complex. Many authors 
found no molecular effect on stopping power. Reynolds et al. [32], and 
Park and Zinmerman [35], found there existed deviation from Bragg's rule. 
Since 1971, Baylor group did several experiments and thus systematical ly 
studied deviation from Bragg's rule due to chemical structure, [6], [7], and 
[8]. They used a particle as the projectile from 300 keV to 2 MeV and 
many gaseous compounds as targets in their experiments. The following 
conclusion have been given: 

1. Physical effects appear to have caused deviation from Bragg's rule. 
The stopping power of H 2 O vapor obtained by Reynolds, et al. [32], 
was found to be an average 11% higher than that of D 2 O ice 
obtained by Wenzel and Waling for proton of 30-600 keV [43]. 

2. Chemical binding effects are more likely to cause departure from 
Bragg's rule for low velocity projectile. 

3. Bourland and Power [7] said Bragg's rule applies to the gaseous 
compounds which contain single and double bonded molecules. Bragg 
rule does not apply to compounds containing triple bonds. For a 
particles of energy 0.3 - 2 Mev, the deviation of Bragg's rule are 
found as much as 12.8%. Especially they indicated that molecular 
hydrogen obeys Bragg's rule. One year later, Power et al. [6] 
wondered about their previous conclusion. They said "this 
observation greatly weakens the assumption that a physical state 
effect is possibly the cause of deviation from Bragg's rule and may 
even imply that the problem is not due to a difference in e(C) 
(stopping cross section) under certain circumstances but rather 
than the atomic stopping cross section e(H) may be considerably 
different than one-half the molecular stopping cross section s (H 2 ) 
as has usually been assumed in the past. In 1974 Lodhi and Power 
[8] gave a more careful conclusion. They said that single bonded 
compounds involving C, H, F and Br have been shown to have molecu- 
lar stopping cross sections that are predictable with errors of a 
very few percent by using vapor deposited solid carbon e(C) along 
with an e(H) that is cornnon to eleven compounds. It appears that 
there exists no unique atomic stopping power for carbon and hydro- 
gen which satisfies Bragg's rule for double bond compounds and that 
due consideration must be given to molecular structure when pre- 


dieting molecular stopping power from atomic stopping power for 
those compounds. 

It appears they corrected their previous conclusion that the double 
bonded compounds have no deviation from Bragg's rule. However, their main 
difficulties lie in the fact that it is very hard to determine the atomic 
stopping power from experimental data. 

Discussion on Deviation from Bragg’s Rule, Due to 
Chemical Structure from Theoretical Model. 

In the above section, we found the Baylor group met difficulties with 
determining the atomic stopping cross section from experiments. Actually, 
it is very hard to obtain atomic state hyirogen. However, from the local 
plasma model it is quite easy to calculate the atomic cross section, since 

* v 

' a atomic electron ground wave functions are employed in obtaining the density 
in local pi asma model . 

Table 11, 12, 13 show the calculated atomic and molecular stopping 
cross section of 0 2 , N 2 , H 2 , respectively. The percentage deviation from 
Bragg's rule are also listed. 


Table 11. Atomic and molecular stopping cross section and deviation from 
Bragg's rule of 0, 


E(keV) 


S( atomic) 
eV + 15"* 5 cm 2 


2 I 15.89 


molecule 


deviation 


0.3 % 2% 


500 

1037 

100000 

9.1 

5.72 

0.1429 

8.92 

5.64 

0.1476 

2 % 

2.4% 

1.1% 


Table 12. Atomic and molecule stopping cross section and deviation from 
Bragg' s rule of N, 


5 (atomic) 

o\i v in-1! 


atom 


2 molecule 


deviation 


100 

200 

18.57 

14.52 

17.79 

1*3. 26 

4.2% 

7.4% 


500 

1037 

100000 

8.53 

5,30 

0.1340 

8.10 

5.20 

0.1319 

5% 

1.9% 

ms 






















































Table 13. Atomic and molecular stopping cross section and deviation from 
Bragg's rule of Hg 


E(keV) 

100 

200 

300 

500 

800 

1037 

2591 

S(2xatomic) 

2xeVxl0“ l5 cm 2 

12.7 

8.13 

6.1 

4.17 

2.89 

2.36 

1.11 

atom 

S(molecule) 

eVxlO“ 15 cm 2 

11.43 

7.53 

5.71 

3.93 

2.75 

2.24 

l.% 

molecule 

deviation 

< 

10% 

7.4% 

6.4% 

. ... 

5.8% 

4.8% 

5.1% 

3.6% 


From Table 11, 12, 13 the following facts are found: 

1. When the projectile's velocity becomes extremely large, the devia- 
tion from Bragg's rule almost vanishes for N 2 , 0 2 gases. When the 
protons energy increased to 100 Mev, the deviation from Bragg's 
rule decreased to almost 1%. This result agrees with Thomson's 
predictions. 

2. When the protons velocity becomes comparable to the atomic electron 
velocity (corresponding to proton energy 40-100 keV) there may 
exist considerable deviations from Bragg's rule. 



3. 


Deviations from Bragg's rule are also found to depend on chemical 

structure of molecule. As is well known in chemistry (b's 
structure is a little bit fuzzy. From bond energy point of view, 

02 is still a double bonded molecule. Meanwhile N 2 is 100% triple 
bonded molecule, H 2 is a single bonded molecule. It is also 

noticed that the maximum deviation for O 2 , N 2 and H 2 , above 100 Kev 
proton energy are about 2,6%, 7.4% and 10%, respectively. Mean- 
while the internuclear distance for Ofe , N 2 , and H 2 are 1.2QA* 

1. 094A and 0. 74A*, respectively . 


It appears, the smaller the internuclear distance the more the devia- 
tion will be. When Gordon-Kim model is employed, the smaller the internu- 
clear distance always means that there is more overlap of electron clouds. 
Thus more overlap causes more deviation from Bragg's rule. In this model, 
the molecular binding effects are also determined by these overlap of 
electron clouds. However, Gordon-Kim model is a very simple model for 
diatomic molecules. The internuclear distance is the only relevant 
parameter but the most important information about the molecular effect is 
contained in this parameter. There is a very strong relationship between 
the bond energy and the distance between the nuclei. The stronger the bond 
energy, the shorter the distance will be. It is interesting to note that 
the single bonded, double bonded and triple bonded carbon molecules have 
inter-nuclear distances equal to 2. 94, 2.52, and 2.24 in Bohr units 
respectively. It means that the triple bonded carbon has more overlap than 
the single bonded carbon. We may thus expect that the triple bond carbon 
will have more deviation from Bragg's rule than the single bonded carbon. 

It is expected from the above statement that for the same compounds the 
triple bonded molecules most likely have more deviation from Bragg's rule 
than single bonded molecules but it does not mean that we agree with the 
statement that the single bonded molecules have no deviation from Bragg's 
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rule. 


Actually from our calculation, it was found that H 2 had considerable 
deviation from Bragg's rule at low projectiles energy. (10% for 100 Kev pro- 
ton energy) due to its small internuclear distance 1.074\* and relatively 
large overlap of electron cloud. It was also noticed that all these devia- 
tions due to molecular effect always decreased the stopping power. 

One of the reasons to understand the above fact is that the binding 
effects always weaken momentum transfer and cause the upper limit of 
momentum transfer to be less than 2mv and thus reduce the stopping power. 
Another reason is due to the shell correction term. As we know the first 
shell correction term is proportion to <T>/mv 2 for bond states <T> ~ 

|<E>| . Thus, more binding effect increased the kinetic energies thus de- 
creasing of stopping power. 


CONCLUSION 


We have established a modified local plasma model which in conjunction 
with the Gordon-Kim model affords a method to calculate molecular stopping 
power even at low projectile energy. By using this model the main conclu- 
sion on deviation from Bragg's rule, summarized by Baylor group's experi- 
ments, can be understood. Some ambiguity on deviation from Bragg's rule of 

is now understood under the present model. The assumption that e (H) = 
j e ( H2 ) , is not correct from overlap point of view (where e is the 
stopping cross section). The only conclusion by using simple Gordon-Kim 
model is that the more percentage overlap caused more deviations. In other 
words, the stronger the binding effect the more deviation from the Bragg's 
rule will be. It appears that for the same compounds, triple bonded 
molecules most likely caused more deviation than the single bonded mole* 
cules. 
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APPENDIX 


THE EQUI PARTITION RULE AND NUCLEAR MOMENTUM RECOILING 

A semi -classical three body collision model has been established to 
estimate the upper limit of total momentun transfer of proton incident on 
hjdrogen atoms. Numerical results revealed that the equipartition rule in 
shell correction deserves more careful study. 

We have mentioned in the section titled "Stopping Number Function" that 
to apply a correction on the first term of shell correction we used Brown's 
results. Specifically the term <T>/mv 2 was used instead of the term 
<T>/|- mv 2 of plasma model. 

In their paper "Stopping Power of Electron Gas and Equi partition Rule" 
Lindhard and Winther [15] mentioned that their result of first term of shell 
correction was in agreement with Walske's result [22]. As is known, 

Walske's shell correction term is just twice that of Brown [21]. This 
result leads one to believe the existence of an equipartition rule in case 
of shell corrections. Lindhard and Winther noticed this, and surmised that 
there was a correspond- ing equipartition rule in the plasma model. This 
implies that the plasma resonance excitation and the close collision each 
had equal contribution to stopping power. It appears that plasma model 
gives exactly the same results as quantum perturbation theory. Lindhard and 
Winther emphasized this fact as a success of local plasma model. Fano also 
mentioned this fact in his paper "Penetration of Proton a Particles and 
Mesons" [44]. He used the results of plasma model to support the existence 
of the equipartition rule in shell corrections. 

The so called equipartition rule has its origin in Bethe* s stopping 
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power formalism. Bethe divided the stopping power of fast charged particles 
into two parts, that due to the distant collision and the close collision. 
There was also found to have considerable contribution to stopping power for 
those particles which move not so fast. Thus besides the logarithm term, 
the so called shell correction terms should also be included in calculating 
the contributions to stopping power. 

Brown had calculated the stopping nunber function and shell correction 
of K shell electron by using the hydrogenic wave functions. The expres- 
sion for stopping number is as follows 


b k 


j Emax^g j Qmax dQ |p (q)|2 
Emin Qmin Q 2 n 


where Q = (P - P') 2 /2m P and P' being the momenta of incident particle 

before and after collision | (Q) | 2 is the form factor. Brown took the 

simple two body collision model, namely particles colliding with a free 

electron to estimate the upper limit of momentum and energy transfer. Thus 

he obtained the maximum momentum transfer A P „ = 2mv. Brown obtained an 

max 

asynptotic expression for stopping power of K electron as & K = 2 In n K + 

1 i 2 • 

2.57861 - — where nj< = — m \r in units of Z e ff Ry is the first term of 
"K 2 K 

shell correction. Meanwhile Walske, took both the upper limit of momentim 
transfer and energy transfer as infinity and obtained the stopping nunber 
functions 


B i ( v n «> * f , EdE f, ^ i F w (Q) i 2 

1 u - a /« 2 ..2 ! A'n n 2 "•! 


Win* 9 j/ n -j 2 w2/4 ni W i 



I 

f 


t 


N 

S 



I •- . 



■ 

it. 


V) 



■j 
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where i ■ K, L, M denotes the different shell, = mv 2 /2 in units of 

Z 2 eff Ry, 0. is the observed ionization potential of ith shell divided by 
"ideal ionization potential." 

Walske also defined the shell correction term C.(0., r\.) as follows 

1 1 1 

(0^ n^) = s. (e.) lnr^ + T. (0^ ) - C i (0 1 , r\.) and he thus obtained 
asymptotic formula B K and both for K and L shell. The correct 
coefficient of l/n K 2 is taken from Khandelwal's paper [25]. 

B K = 2 In n K + 2.57861 - 2njj 1 - (£i) n |< - 2 
B^ = 8 In n^ + 25.5766 - 2n^ (to order n^ 1 ) 

Notice that the first shell correction terms are 2H- 1 or 20 - 1 . The K- 
shell term 2 n^ 1 thus is twice that of Brown's term n- 1 . 

In other words, shell correction can also be divided into two high and 
low momentun transfer parts. Such parts each contributes equally to 
stopping power. Indeed, Walske explicitly divided shell correction into two 
parts. 

C (0, n) ■ Ci (0, n) + C 2 ( 0 , n ) 

where Ci and C 2 are low and high momentum transfer parts respectively. 
Furthermore, he also showed that Ci and C 2 are equal to the order of nf 1 
both for K and L shells and to order l/n K 2 for K-shell. Walske' s 
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results are only limited to K and L shell electrons of hyirogenic atoms. 
But Fano made an assumption that the conclusion Ci s C 2 can be generalized 
to any shell of any atom. This is so called the equipartition rule of shell 
correction. All these conclusions are based on Walske's assumption that the 
upper limits of the momentun transfer and the energy transfer are infinity. 
The question is whether the assumption is true or not. 

Brown took the upper limit of momentun transfer as 2mv and this is 
based on the assumption that the electron initially is free. Two factors 
are neglected in this assunption. First Brown neglected the binding effect 
of the electron. Second, he also neglected the nuclear momentum recoiling. 
It is obvious that if the nuclear motion is involved then due to its huge 
mass, the total momentun transfer may be greatly increased. But by how 
much? Is taking Infinity a good approximation for the upper limit of mo- 
mentum transfer? It appears the answer should be dependent on the velocity 
of projectile. For instance, if the projectile moves extremely slow then 
the nucleus may obtain sufficiently large momentun transfer. Otherwise due 
to the short interaction time, nuclear momentun could be small and the upper 
limit of total momentum transfer will not differ too much from 2mv. 

In this appendix, we estimate the upper limit of momentun transfer by 
the projectile to a hydrogen atom by establishing a semi-classical three 
body collision model. For protons energy over 50 Kev a fitting formula from 
numerical results was obtained: (See for the general results later in this 

appendix) 


AP = 2 m V 5 
max e p 
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where 5 » | (1 + 8.9 (^)S) (1 + J 1- (J£)2) 

V P " V P 

» 

where V p is the velocity of proton and V Q is the Bohr velocity. Nuneri 
cal results of £ of various proton energies are listed in table 14. 


Table 14. 


E(keV) 

50 

75 

100 

200 

400 

1600 

£ 

2.20 

1.43 

1.19 

1.01 

0.991 

0.996 


From the above formula and the table we see that for high energy pro- 
jectiles Brown's assumption is correct. The values of £ slightly less 
than unity is due to binding effect but for low energy proton the factor of 
nuclear momentum play a more important rule. There is a considerable 
correction to Brown's results. But even at 50 Kev protons energy the 
numerical result of £ is only 2.2, st i, 11 quite different from Walske's 
assunption of infinity. It is true that when protons energy becomes 
smaller, then the correction factor £ is expected to increase rapidly, but 
then the expansion of stopping number function B^ should involve more 
terms than just the term n^” 1 . However Fano' s assumption that Cj - C 2 or 
equipartition rule of shell correction, deserves a careful study. 

Nuclear Momentun Recoiling 

In this section first we shall review the conservation laws and the 
resulting physical quantities when the collision is assumed to be taking 





place between two bodies only. Later, three body collision problems are 
handled numerically in a semi-classical manner. Nuclear momentum recoiling 
effect is estimated in the relatively low energy region of Bethe stopping 
power. 

Semi-Cl assical Formulation of Three Body Problem for Relatively 

Low Energy Stopping Power 

Estimating of the Quantities Q(min) and Q(max) Two Body Collision . 

Let us first consider the energy regions such that we can always ne- 
glect the nuclear recoiling energy, i ,e. we can express energy conservation 
law as 


(E„ - E 0 ) - h 2 (K 2 -K' 2 )/2M p 

where E n and E Q are the eigen energy values of final state of electron 
and initial state of electron, hK and hK ' are the initial and final momen- 
tum of projectile. 

We know that in lab. coordinates the momentun changes of projectile are 
much smaller than the momentum of projectile itself, i.e. aK<<K or 0<<1 so, 
we have (E n -E Q ) = h2(K2-|<2 ' )/2M s h 2 K(K"K ' ) /Mp =h V p (K*K ' ) 

since 0 is small 

Q2=K 24 K' 2 -2KK' cos(0 )=(K-K') 2 + (K e ) 2 ' 

q *V (E n- E o /h V 2 + V 2 



Thus Q(min)=(E n -E 0 )/hVp 


This Q(min) estimation is justified unless one deals with the extremely low 
energy case for hydrogen atom targets when proton's energy is lower than 10 
kev (Andersen and Ziegler) for which we need to consider nuclear recoiling 
energy term. Now let us estimate Q(max) term very carefully. This upper 
limit estimation will involve a correction to Bethe's theory. 

We use Landau's treatment of Bethe' s theory for estimating the upper 
limit of momentum transfer. Essentially it neglects the nuclear recoiling 
displacement and momentum. Hence it becomes simply a two body collision. 

In this model, the projectile collides with a free electron. 

The projectile's initial momentum is hK 
and the final momentum is hK'. All the 
momentun changes of projectile are 



transferred to the electron which is free in its initial state henc$ we 
writedown the momentum and energy conservation laws as. 


K' 2 

AE 


= K 2 +aK 2 -2KAKcos<j> 

k 2 -k* 2 . q 2 

- ar«MfUP 

2 m m 

p 


AK * Q 


where Q is the momentum of electron (in our case projectile momentum loss 
is entirely transferred to the electron) m is the mass of electron. From 
above conservation laws, we have 
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2KAKcos4> -AK 2 _ Q 2 


2M 2m 

2KQCOS9 a (1 + 1_) q2 * 

M„ m m 

P P 


or Qf max) = 2 mV p 

It is a very simple model, but it works at high energy. The main physical 
reason is that the projectile moves too fast to cause any significant 
displacement of the nucleus. At relatively low energy due to long 
interaction time, we need consider the nuclear recoiling momentum to correct 
the upper limit i.e. Q(max) 


Semi-Classical Three Body Formulation 
A model essentially based on classical three body problem was estab- 
lished. Consider a proton which collides with a hydrogen atom. The 
electron of the atom was initially assuned to be in its ground state. Some- 
what semi-quantim conditions are involved in our model, besides usual 
classical colunb interaction. Now let us see how the nuclear recoiling 
momentim influences the upper limit. 



where aP q and aP r are the difference of the electron momentim and the 
nuclear momentim before and after collision. Now the projectiles momentim 


<C - 3l. 
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OkS'-itiViU | *'i ' ' JT: 

OF POOR QUALITY 

not only can be transferred to the electron but also can be transferred to 

the nucleus. Furthermore as stated before we assune that the electron is in 

the ground state initially. Now let us set the initial momentun of nucleus 

to be zero i.e.» AP =p n . Thus we have AK=AP +P„ 

n r n e n 


. (aP ) 2 + 2aP + (P n ) 2 

LK - |PJ /_? _i-2 5— • |P 1(1 + n) 

1 


where 1 + n = 


(AP e ) 2 + 2AP e -P n + (P n ) 2 


P 2 

e 


( 20 ) 


We again write down the conservation laws; 


K‘ 2 = K 2 + AK 2 - 2KaK co s4> 

K 2 K' 2 Pe 2 Pn 2 

1 JL*. + e = 2 + _ + _ + V 

) 2M 2M 2m 2M n 1 

P P n 

AK = |P | (1 + n) (21c) 

^ e 

For maximim momentun transfer cos(<|») = 1, as before. From (21a) we have 

K 2 -K 2 '/2M =(2KaK-aK 2 )/2M 
P P 

from (21a), (21b), (21c) 
we have 

p 2 p 2 p 2 

V„|P.| (1 +n) - (1 +n) 2 =_L. + JL- (E 0 - V ) 

P 2M„ 2m 2M„ n 

p n 


(21a) 

(21b) 
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hence, 


ORIGINAL PM® 

OF POOR QUALIFY 


v pl p el (1 + "> + < E o - V„) 




where a ■ + 1__ A 1 P n 

M m M P 2 

P n e 


or p e 2 ** 



V 1 * 1 ) / (V -E )a 

P e - 1 - 2 - 2 - 


( 22 ) 


V (1+n ) 2 j (V -E )a 

'■ 4I W ■ '.a*.) - a J i. .‘W ) 

2 p 


(23) 


For high energy case V p »V Q and thus there is no nuclear recoiling 


cH-l/m, n*0 (V n -E 0 )o/ iv p (l+n)2> 0 


hence = 2V p m 

Therefore, for relatively low energy case, we have two corrections; one is 
from nuclear recoiling momentun, and another correction is from the initial 
energy of electron. For most cases the first correction is more important. 


OF FOOtt QbALi Vi 

from the formul a for AK max , we can see AK max depends on 1+n or depends 
on nuclear recoiling momentum P n which extends the upper limit and will 
cause an additional stopping power. 


In formula 1+ n 


AP a 2 + 2 aP P n P 2 
e e n + n 




(24) 


All these AP gt p , and P Q are numerical results of three body problem. 
Now let us establish fundamental equations. 



we denote projectile by I, electron by 
2, and nucleus by 3. 

For simplifying the calculations we established a two dimension model. 

It is not a bad approach because we are interested in nuclear recoiling 
effects in which the main contribution is due to the colutib interactions 
between projectile and nucleus. This is partly due to the symmetry and 
partly due to our interest in obtaining average values. 

Equations of motion for projectiles are as follows: 


M 


M V[ = - 


e 2 Hi -*2 ) + e 2 (Y t -Y 3 ) 

3 + 


12 


13 


« Xj = 


£lhihl ftthjM 


12 


13 
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i 


Yi 


t=to 


= b 


ORIGINAL PAQS IS' 
OF POOR QUALITY 


Xi 


t=to 


4ao 


Xi 


= Vp 

t=to 


Where M is the mass of projectile which in our case is a proton. We took 

the integration time as projectile moved from -4ao to 4ao for the standard 
gas state given the average distances between molecules are about the order 
of 10A and b is the impact parameter. 

For the electron we have the following equations: 


e 2 (Yi -v 2 ) e 2 {Y 2 -V 3 ) 

' ' ' ~ ~ r- r — ^ 


m Y 2 = 

r 3 

* 

r 3 


t=to 

12 

23 



.. e 2 (X 3 -X 2 ) 

mX2 - - 

r 3 

e 2 (X 2 -Xi ) 
r 3 

X 2 

= -COS^) 

t=to 

23 

12 



where <f> is the initial 

phase of electron orbit 

Equations of motion of 

nucleus are: 



e 2 (Y 2 -Y 3 } 

e 2 (Yi -Y 3 ) 

y 3 

= 0 

M y 

n 3 r 3 

r 3 

0 

t=to 

23 

13 



.. e 2 (X 3 -X 2 > 

e 2 (X 3 -Xi) 

X, 

f = 0 

«n X 3 ’ r3 ' 

r 3 

3 I 

jt=to 


23 13 


Y 2 = Vo si fy 
t~to 
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As is well known nunerically high order differential equations can be reduc- 
ed to a first order equation such as: 

Y (m) *f(x,Y,Y',y',....y (m “ 1) ) 

Let Y*Yi, Y'*Y 2 , Y"=Y 3 ,....Y^ m " 1) *Y m then this will reduce to 

r Vv 2 

Y 2 '=y3 

4 

, Y , =f(x,Y,,Y 2 ,...Y ) 
v m 12 m 

So our six second order differential equations are reduced to twelve 
first order equations. The electron is assuned initially to be in the first 
Bohr orbit. Since we took laboratory coordinates the nucleus is initially 
at rest in the position of origin. We emphasize here, the difference from 
the usual treatment. Here we took the impact parameter b as the vertical 
distance from the projectile to the nuclei's initial position and not the 
vertical distance from the projectile to the electrons initial position, as 
is usually done. To apply classical mechanics to the microscale system, the 
most serious difficulty is that the electron can eventually drop into 
positive ions colunb potential well. It also caused nunerical difficulties 
in practice because we need infinitesimal steps to keep acceptable 
accuracy. 

To prevent this difficulty it looks as though we need to introduce 
somewhat semi-quantim condition. In 1951 David Bohm suggested an interrup- 
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t Ion of quantun theory in terms of "hidden variables." [45] Bohm proved 
that quantum mechanics can be explained as some modification to classical 
mechanics. Schordinger' s equation can be expressed as: 


dt 2 2m R 

Where R is the real part of wave function or R 2 = 1 (r) | 2 

So the equation of motion for quantum mechanics can be expressed as us- 
ing classical potential plus a quantim mechanical potential which is correc- 
tion to classical theory. 

Following Bohm we got some hint that this semi-quantun mechanics treat- 
ment prevents the electron from dropping into the nucleus. If given an 
additional semi-quantun potential then we will have an additional force to 
balance the usual colunb force. 

For r<<a 0 this force will be greater than columb force. It will pre- 
vent electron dropping into the nucleus. For r>ao, this force will 
vanish. In some sense this additional potential gives explanation of first 
Bohr orbit. The above discussion is for the bound state. For the 
scattering state the electron can approach the nucleus. Thus, we will 
establish a potential somewhat like Fermi distribution which depends on the 
energy of the electron 


ed 


(e- 10 +l) 
rf (e lCE +1) 


(25) 


Figure 2 gives the graphics of this potential for £ >1 Ved~0, and for E<-1 
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V or we /nay say for scattering state our model Is a classical 

ea 2 mR 2 

model but for bound /State we have an additional term somewhat semi-quantum 
term. 

We use the expression (25) rather than a potential like 
v . h 2 E<0 

ed 2mR 2 


V fid = 0 E>0 

merely due to the reason that we need a continuous potential to obey the 
conservation law accurately. We have some freedom to choose e- 10 or e- 12 
or e- 8 in the terms of V e(j . The bigger the number we choose in 
exponential terms the sharper the curve obtained. But as mentioned too 
sharp a value will cause energy or momentum conservation problems. 

We calculated these recoiling momenta and averaged them over phase of 
the electron orbit and then took the average over impact parameter b. Here 
as we mentioned earlier, b is the vertical distance from projectile to the 
initial position of the nucleus. From (23) we have the maximum momentum 
transfer 


AK = V iltlll 2 - 
max p 



< v n*o>« j s 

^ V (1+n)2 


? 2m v p 


In the stopping number formula instead of 2m Vp we have 2rn Vp? hence we 
finally obtained the new stopping n umber formula 
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In this two dimension model, there are two parameters, impact parameter b 

and initial phase. Nunerically we calculated P n , P g , Pp, (1+n), and 

finally took the average over 4> and b. Table 15 shows some examples of 

P , P , P as functions of initial phase. Table 16 shows some examples of 
n 6 p 

the average of <p as a function of impact of parameter b. The correction 
coefficient of the stopping nunber involves twd factors. One is an addi- 
tional contribution due to nuclear momentim recoiling. As b is very small 
this factor is very important, as b becomes larger this effect vanishes. 
Another factor is mainly due to the velocity of the electron in the initial 
state. This factor makes negative contribution on correction (make it less 
than unity) as V p increases this factor becomes negligible. After an 
average over <f> and an average over b, we obtained table 17 as a function 
of the projectiles energy. 


Table 15. Momentuns after collision as a function of the initial phase. 


Initial Phase 

Momentums after collision 

♦ 

Pp(me-Bohr/sec) 

P e (m e - Bohr/sec) 

P n (m e -Bohr/sec) 

0 

3672 

1.39 

9.28 

1.10 

3672 

1.36 

11.47 

2.04 

3672 

1.30 

11.27 

2.98 

3667 

2.80 

2.86 

4.08 

3671 

1.83 

12.31 

5.03 

3667 

4.31 

7.65 

5.97 

3672 

2.35 

9. 67 


(a) b ».i (a 0 ) E = 100 kev 


Initial Phase 

Momentuns after c 

:ollision 

' 

4> 

Pp(m e -Bohr/sec) 

P e (m e -Bohr/sec) 

P n (m e -Bohr/sec) 

0 

3672 

1.47 

1.99 

1.10 

3672 

1.37 

3.63 

2.04 

3672 

1.21 

3.50 

2.98 

3671 

1.73 

2.88 

4.08 

3671 

2.15 

0.84 

4.87 

3668 

4.72 

0.95 

5.97 

3672 

1.25 

1.09 
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As we expected when projectiles energy Is small that is when the pro- 
jectile moves slow, it causes significant nuclear momentum recoiling effect. 
(Meanwhile the initial velocity effect is covered by nuclear recoiling 
effect. As projectile moves fast it can not cause significant nuclear re- 
coiling momentun. We can see that the Initial velocity effect however is 
less important when the projectile moves very fast both effects vanish and 
there is no correction and we obtain the Bethe's formula. 

An approximate useful forinula obtained by fitting the values as energy 
greater than 50 kev, is: 


5 


1 

2 


(1 + 8.9 (— ) 5 ) (1 
V P 


+ Jl - (— ) 2 ) 

' V P 


Where V Q is the Bohr velocity and V p is the velocity of projectile. 
Numerical results relative correction of stopping number are also calcu 
ed. 

L/Lq as a function of projectiles, energy are listed below 


E(keV) 

50 

75 

100 

200 

400 

1600 

L/Lq 

1.398 

1.15 

1.065 

1.0035 

0.998 

0.9993 


From the table we can see that for the hydrogen atom target and for the 
proton as a projectile of 50 keV nuclear recoiling momentim caused a consid- 
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erable correction (40%) on stopping power for 100 kev about 6-7% correction 
on stopping above 200 kev the correction of stopping power can be 
negligible. 


Argument of Extending to General Material. 

Up to now we have calculated the recoiling momentun of hyderogen atom 
target. Now we will give an argument that this result can be roughly 
extended to a general case of any other atoms. Let us consider proton 
projectile passing through material composed of atoms of charge Ze. Then 
due to the colunb interaction between proton and nucleus the momentum trans 
fer is: 


iP n * C y *> « 

where f (t) is proportional to Zt and is independent of the nuclear 
mass. Now only for estimating we suppose that all the momentim obtained by 
nuclei is transferred to the electrons. Thus Z electrons share their 
additional momentim. '-fence on the average each electron obtained momentim 
independent of Zt and nuclear mass. 

Then the additional stopping power is proportional to Zt but the 
stopping nunber is independent of Zt and nuclear mass. We can roughly say 
that the ratio of "additional stopping nimber" to stopping number of any 
atom are the same as for the hydrogen atom. 


